clear
close all
clc

cd ..\
cd ..\
cd Stored_file

load IRF_demo

cyshare=DEBTIRFn;
cysharefb=[Ecyfb Ecyfb_shock Ecyfb Ecyfb Ecyfb Ecyfb Ecyfb Ecyfb];
CYsharefb=[ECYfb ECYfb_shock ECYfb ECYfb ECYfb ECYfb ECYfb ECYfb];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%% FIGURE D1 %%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

plot([1:1:8],cyshare,'k','LineWidth',1); hold on;
scatter([1:1:8],cyshare,'filled','MarkerFaceColor',[0 0 0]); hold on;
plot([1:1:8],CYsharefb(1).*ones(1,8),'k:','LineWidth',0.1); hold on;
scatter([1:1:8],CYsharefb,'filled','MarkerFaceColor',[0.5 0.5 0.5]); hold on;
plot([1:1:8],cyshare(1).*ones(1,8),'k:','LineWidth',0.1); hold on;
plot([1:1:8],CYsharefb,'color',[0 0 0]+0.7,'LineWidth',1);
set(gca,'FontSize',09);
xlabel('Time')
ylabel('E[c],E[c$$^{\ast}]$$','Interpreter','Latex','FontSize',10)
box on

